Evidence for "fragile" glass-forming behavior in the relaxation of Coulomb frustrated 

three-dimensional systems. 
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We show by means of a Monte Carlo simulation study that three-dimensional models with long- 
range frustration display the generic phenomena seen in fragile glassforming liquids. Due to their 
properties (absence of quenched disorder, physical motivation in terms of structural frustration, and 
tunable fragility), these systems appear as promising minimal theoretical models for describing the 
glass transition of supercooled liquids. 
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The viscous slowing down of supercooled liquids that 
leads to glass formation when the temperature is suffi- 
ciently lowered, remains a challenging problem of liquid- 
state and condensed-matter physics. The phenomenon 
seems dramatic enough (a change by 15 orders of mag- 
nitude of the viscosity rj and the primary (a) relaxation 
time Tq, for a mere decrease of the temperature T by a 
factor of 2) and general enough to call for a universal ex- 
planation. Yet, the search for an underlying universality 
is thwarted by the absence of observed criticality, e.g., 
by the absence of a diverging, or rapidly growing, static 
susceptibility or correlation length j^]. 

If one accepts the premise that the distinctive fea- 
tures of most glassforming liquids, namely the rapid 
super- Arrhenius increase of rj and with decreasing 
T and the non exponential character of the relaxation 
functions, result from collective or cooperative behav- 
ior, progress could be made by identifying the key phys- 
ical ingredients and formulating a minimal theoretical 
model based on these ingredients, in much the same 
way as spin glasses have been described in terms of 
quenched disorder and frustration and modelled by the 
Edwards- Anderson Hamiltonian. One such attempt has 
been recently madeQ, based on the concept of structural 
frustration!^. Such frustration arises because a locally 
preferred structure, i.e., an arrangement of neighboring 
molecules that minimizes some local free energy, cannot 
tile the whole space, because of the global constraint!^. 
The canonical example is a 3-dimensional system of 
spherical particles whose locally preferred structure, an 
icosahedral cluster of 13 particles, cannot form a pure 
crystal because of the 5-fold rotational symmetry of the 
icosahedronQ. In the frustration-limited domain the- 
ory (FLDT) of supercooled liquids|^, structural frustra- 
tion is implemented as the competition between effec- 
tive interactions acting on very different length scales: 
a short-range ordering term that favors the extension of 
the locally preferred structure and a weak, but long-range 
(l/r) frustrating term that generates a strain free energy 
growing super-extensively with the linear size of the or- 



dered region. These ingredients are most simply incor- 
porated in three-dimensional Coulomb (l/r) frustrated 
lattice models with the following Hamiltonian, 
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where J and Q are both positive and denote the strength 
of the ordering and frustrating interactions, respectively, 
< ij > is a sum over distinct pairs of nearest neighbors 
and \rij \ is the distance between sites i and j (the lattice 
spacing is taken as unit length); the tunable ratio Q/J 
characterizes the frustration strength. In such models, 
glassiness has to be self-induced and does not result from 
quenched randomness or from specific kinetic constraints. 

In this letter, we present an extensive Monte Carlo sim- 
ulation study of Coulomb frustrated models on a cubic 
lattice. We focus on the relaxation to equilibrium and we 
do not consider at this stage out-of-equilibrium dynam- 
ics and aging phenomena in a non-equilibrated, glassy 
state. We show that: (i) the relaxation displays the main 
features characteristic of the slowing down of "fragile" 
glassforming liquids (non exponential decay of the relax- 
ation function, super- Arrhenius activated T— dependence 
of Tq); (ii) the "fragility" of the system, e.g., its degree of 
departure from Arrhenius behavior, increases as the frus- 
tration decreases, thereby providing a description that 
spans the whole range from strong to very fragile glass- 
former via tuning of the relative amplitude of the frus- 
trating interaction; (iii) the slowing down occurs without 
the rapid growth of a correlation length and does not 
follow the standard critical slowing down pattern; (iv) 
the results, including the existence of a crossover tem- 
perature in the vicinity of the critical point of the un- 
frustrated system, confirm the scenario proposed by the 
FLDT. 

In using Eq. (0), we have chosen spin variables Si 
that can take several orientations; most of the results 
have been obtained for the 5-state clock model (Sj = 
{cos{2Tr9i/q),sm{2TT9i/q)), with 6i the orientation of the 
planar spin and q = 5) and the Ising model {q — 2), but 
we have also studied the g = 11 state clock model. 

In the absence of frustration (Q = 0), the models have 
a critical point at a temperature T^, below which they are 
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FIG. 1: Spin-spin correlation function C{t) versus ln(t) for 
several temperatures: 5-state clock model for Q = 0.1 (a) 
and Q = 0.00625 (b), 11-state clock model for Q = 0.05 (c) 
and Ising model for Q = 0.001 (d). In all cases, curves from 
left to right are for decreasing T's. 



ferromagnetic. The long-range, coulombic interaction re- 
quires that the total magnetization of the system be zero 
in order to ensure a proper thermodynamic limit. There- 
fore, long-range ferromagnetic order is prohibited at all 
T's for any nonzero value of the frustration parameter 
Q/J; in the FLDT, this mimics the effect of structural 
frustration that forbids any long-range order based on 
the periodic repetition of the locally preferred structure. 
The models with discrete orientations studied here can 
still form ordered phases characterized by modulated pat- 
terns (e.g., lamellar phases) Q, but as shown recently by 
Monte Carlo simulation of the Ising model ||], the transi- 
tion between paramagnetic and modulated phases is first- 
order. Following BrazovskiiQ, this result can be inter- 
preted on the basis of the self-consistent Hartree approx- 
imation which predicts the occurence of a fluctuation- 
induced first-order transition, a transition with no nearby 
low-T spinodal. All the simulations discussed below are 
made in the disordered (paramagnetic) phase. Since the 
transition to the modulated phases at TooiQ) first- 
order, one could in principle supercool the paramagnetic 
phase to lower temperatures; however, we have found 
that the lattice sizes achievable in practice are too small 
to allow for a proper supercooling below Too- 

To study the relaxation of the frustrated systems, we 
have considered the dynamics associated with the Monte 
Carlo algorithm (simple or slightly modified Metropo- 
lis rule0]), time being the number of sweeps per spin. 
In all cases, the total magnetization is constrained to 
be zero, as it should be in the thermodynamic limit. 
The dynamical quantity that wc have monitored is the 
spin-spin self correlation function, C{t) — {l/N)J2i < 
Si{t')Si{t' + t) >, where the bracket denotes the ther- 
mal average and N is the total number of lattice sites. 
We stress that C{t) is an equilibrium correlation func- 
tion that is computed once the system has reached equi- 
librium. In practice, the thermal average has been per- 




FIG. 2: C{t) versus ln(t) for the 5— state clock model with 
Q = 0.00625. Left curve: T = 3 > T°; the dashed line, virtu- 
ally indistinguishable from the data, is an exponential. Right 
curve: T — 1.7 < T^; the dashed curve is an exponential and 
the dotted curve is a stretched exponential with /3 = 0.4. 



formed as an average over 20 different initial times t', 
all chosen longer than the relaxation time. Most sim- 
ulations have been done on cubic lattices of size 16"^ to 
22^^ with periodic boundary conditions, and the Coulomb 
interaction has been handled via Ewald sums. In one 
special case (see below), we have performed a finite-size 
study with lattices of linear size L = 5, 10, 15 and 20. We 
have covered two orders of magnitude in the frustration 
parameter Q/J, always in the weak-frustration regime 
(Q/J << 1). In what follows, we take J —\. 

The evolution of C{t) with T is illustrated in Fig. |l| 
for the 5-state clock (Q = 0.1 and Q = 0.00625), the 
11-state clock (Q = 0.05), and the Ising (Q = 0.001) 
models. At the highest frustration, illustrated by Fig. |^a, 
the decay of C{t) appears to proceed in a single step at 
all T's. For lower frustrations, a 2— step decay develops 
as T is lowered (see Fig. |l|b and d). At high T's, typically 
above the critical point of the coresponding unfrustrated 
model (T° ~ 2.1 and r° ~ 4.51 for the 5-state clock 
model and the Ising model, respectively), the whole time 
dependence of C{t) is well fitted by a simple exponential, 
exp(— t/To(T)); this is shown in Fig. ^ but is true for all 
models and all frustrations. At low T's, below T°, it is 
impossible to describe the entire decay of C{t) by a single 
exponential; as illustrated in Fig. ||, the emerging second 
step of the relaxation can be described by a stretched ex- 
ponential, eyip{—{t/TKWw{T))^'^'^^), while the first step 
can still be fitted by a simple exponential, exp(— i/ro(T)). 
This non-exponential behavior and emergence of a two- 
step decay, both of which become more marked as T is 
decreased, are typical of many actual glassforming sys- 
tems, especially the fragile supercooled liquids[0. One 
may however notice that the timescale separation be- 
tween the two relaxation steps that can be achieved in the 
simulations is not sufficient to observe the development 
of a true plateau at intermediate times. 

To study in more detail the characteristics of the slow- 
ing down of the relaxations as one lowers T, a phe- 
nomenon already visible in Fig. nl we have monitored 
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FIG. 3: Arrhenius plots of ro (lower curve) and (upper 
curve) for the 5— state clock model with Q = 0.00625 (a) and 
g = 0.05 (b). 



as a function of T the relaxation times characteristic 
of the short- and long-time behaviors, tq{T) and Ta(T) 
respectively. The former is obtained from the fit to a 
simple exponential (see above) , but an equivalent behav- 
ior is obtained from the condition C(to) = 0.9. The 
latter is defined, as often for spin systems, by the con- 
dition C(tq) = 0.1; in this case indeed, there was too 
much uncertainty in using the three-parameter stretched 
exponential fit that also requires the introduction of a 
lower cut-off time. The results are illustrated in Fig. ^ 
for the 5-state clock model. We have plotted the log- 
arithm of To and versus 1/T. One can see that 
over the range of T studied, the dependence of to(T) 
is essentially Arrhenius-like, i.e., linear on the diagram, 
To(T) ~ To,oo exp(i?oo/T), whereas the dependence of 
Tq(T) shows a marked deviation from Arrhenius behav- 
ior below some crossover temperature in the vicinity of 
T°. Note that at high T's (above T°), C{t) is essentially 
exponential and the two times to{T) and Ta(T) differ by 
a trivial constant, so that in this T range the Arrhenius 
dependences of tq and are the same. At lower T's, 
To is characteristic of a "secondary" or "precursor" re- 
laxation; one expects it to be associated with a local, 
weakly or non cooperative, mechanism of relaxation. On 
the other hand, the super- Arrhenius rise of Tq versus 1/T 
is indicative of cooperative behavior. 

The crossover from Arrhenius to super-Arrhenius be- 
havior of Ta is shown in Fig. ^ for the 5— state clock 
model for 5 different frustrations from Q = 0.1 down 
to Q = 0.00625. A similar trend is observed for all 
other models studied (for the Ising model, see Ref. (?[ ||)). 
Several points are worth noting: first, such a super- 
Arrhenius behavior is typical of the viscous slowing down 
of fragile glassforming liquids, the more fragile a liquid 



the more pronounced the super-Arrhenius character ||lj; 
secondly, the crossover occurs in the vicinity of the criti- 
cal point of the unfrustrated system; thirdly, the depar- 
ture from Arrhenius behavior, i.e., the "fragility", be- 
comes more marked as frustration decreases: see Fig. ||. 
To our knowledge, such a "tunable" fragility that pro- 
vides a description over the whole spectrum of glassform- 
ing behavior, from strong to very fragile, is shown by no 
other microscopic model. 
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FIG. 4: Arrhenius plot of Ta{T) for the 5— state clock model; 
from left to right: Q = 0.00625, Q = 0.0125, Q = 0.25, 
Q = 0.05, Q = 0.1. The full lines are the best fits to Eqs. (2,3). 

The competition of interactions acting on very differ- 
ent length scales can thus generate fragile glassforming 
behavior similar to that observed in actual supercooled 
liquids. In the FLDT|^, this behavior is attributed to 
the proximity of the critical point T^ of the unfrustrated 
system, a critical point that is avoided in the presence 
of frustration. Avoided critical behavior||^, ^ results in 
the appearance at T's below T° of domains whose size 
is limited by the frustration. These domains are respon- 
sible for the non-exponential character of the relaxation 
and the super-Arrhenius, activated T— dependence of Tq. 
Scaling arguments predict for instance that the effective 
activation free energy for the a relaxation behaves as j^] 



E^{T) = E^, T>T*, 



- i;^ + ST* ( 1 - _ 



8/3 



(2) 

T < T*, (3) 



where Ea{T) = Tin 



AT) 



with Too a characteristic 



high-T time, and T* is close to T°. The existence of 
a crossover in the vicinity of T^ is indeed supported 
by the Monte Carlo results. More quantitatively, we 
have compared the above prediction with the simulation 
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data. In doing so, one must remember that the range of 
Tq's spanned by our data is hmited, as in most simula- 
tions of glassforming systems, and is far less than the 15 
decades observed experimentally in supercooled liquids. 
The high-T parameters Too and E^o are independently 
obtained by an Arrhenius fit to to{T) (see above, where 
T(x,/to,oo = — In(O.l) due to the definition of Tq, and tq), 
and one is left with two adjustable parameters, B and 
T*, for each frustration (and each model). The resulting 
curves are shown in Fig. ^ for the 5— state clock model. 
For all models the agreement between theory and simula- 
tion is good. As predicted by the theory, T* is found close 
to TO (e.g, T*/T° ~ 1.14 - 1.20 for the 5-state clock and 
1.04 — 1.12 for the Ising model) and the super- Arrhenius 
parameter B varies roughly as Q~^. (As already stressed, 
small frustration Q implies large super-Arrhenius char- 
acter and large fragility.) 
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FIG. 5: Arrhenius plot of Tq for the 5— state clock model and 
for lattice sizes L = 5,10,15,20 (from right to left): Q = 
0.00625 (a) and Q = (b). 



We have thus found that our simulation data are 
compatible with an activated-type expression, Tc(T) oc 
exp(£'„(r)/T) with E„(T) given by the FLDT expres- 
sion. (As discussed in |g, |l^ , the data are also compatible 
with Ea{T) described by the entropic droplet /random- 
first-order transition approach (EDA)[|0[ y, which of- 
fers a unique opportunity for comparing these two al- 
ternative theories in more detail.) On the other hand, 
we have checked that they are poorly described by a 
power-law dependence, T^iT) oc \T — Tc\~'^ , except of 
course when Q = 00. In both the FLDT and the 
EDA, the slowing down of the relaxation, as well as the 
non exponential character is attributed to genuinely non- 
perturbative effects, such as heterogeneities (domains or 
droplets) and cooperative, activated processes. This is 
quite different than standard critical slowing down, as 
could occur if driven by the approach of a spinodal. In 
particular, the length scale associated with the slowing 
down is predicted to grow slowly with decreasing T, at 
a much slower rate than the relaxation time itself. To 
check this, we have performed a finite-size study of the 
5-state clock model for Q = 0.00625. As seen in Fig. ||a, 
Tq (T) is indeed unchanged when one decreases the linear 
size of the system from 20 to 10; on the other hand, it 
is strongly altered when passing to L = 5. This seems 
to confirm that any putative length scale associated with 
the strong slowing down is rather modest (of the order 
of 5 or so) over the whole T— range studied. For com- 
parison, the same finite-size study for the unfrustrated 
{Q = 0) model, showing the usual pattern associated 
with critical slowing down, is displayed in Fig. 

In summary, due to their properties (absence of 
quenched disorder, physical motivation in terms of struc- 
tural frustration, 3— dimensional, finite-range character, 
tunable fragility), whose combination is found in no other 
models ll^, the Coulomb frustrated systems appear as 
promising minimal theoretical models for describing the 
slowing down of relaxations in supercooled liquids. 

We thank D. Kivelson and L. Berthier for many stim- 
ulating discussions. 
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